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Abstract. It is well known that Yukawa potentials permit bound states in the Schrodinger 
equation only if the ratio of exchanged mass to bound mass is below a critical multiple of 
the coupling constant. However arguments suggested by the Darwin term imply a more 
complex situation. By studying numerically the Dirac equation with a Yukawa potential 
we investigate this "amplification" effect. 
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I. INTRODUCTION 

In recent years neutrino oscillation experiments have all but proved that neutrino mass states 
exist [1]. These are linear combinations of neutrino flavor states. The neutrinos arc subject only to 
weak interactions (neglecting gravitation) . More precisely, they participate only in weak vertices (loop 
diagrams in which a neutrino couples to a charged lepton virtual pair, W^^ would of course allow for 
higher order electromagnetic interactions via the charged virtual particles). 

Since neutrinos have mass and couple via exchanged intermediate mass bosons to leptons one may 
legitimately ask if a neutrino-lepton bound state could exist. Exchange of , Z° means the treatment 
of Yukawa potentials and the question of when bound states exist for Yukawa potentials [2] . Such weak 
bound states would provide very interesting effects in superconductivity as they would constitute a 
"bosonic atom" which remind us of one of the multiple roles played by Cooper pairs in standard 
theory. 

The question of the existence, or not, of a bound state for a given potential seems a simple 
theoretical question with a straightforward method for flnding the answer. Given a potential and the 
corresponding reduced mass wave equation one solves it, normally numerically, and sees if bound 
states exist. Bound states are characterized by normalized (localized) solutions and a discrete energy 
spectrum below the free particle threshold E < m ^ iS^R = E — m < 0, where E is the relativistic 
energy and iS^R the non-relativistic energy of the particle with reduced mass m. Either the existence 
of normed states or of a discrete spectrum suffices to identify a bound state regime. Of course, nothing 
is quite so simple. For example there are limits for the validity of the use of potentials. Furthermore 
coupling constants have the "annoying" tendency to run and hence, are all but constants. 

Even within the realm of non-relativistic quantum mechanics potential theory, we have surprising 
subtleties. A one-dimensional square well always yields a bound state. A three-dimensional spherical 
well must be sufficiently "deep" or extended to allow a bound state solution [3]. Relativistic wave 
equations open up numerous conundrums. Can a bound state energy exist below the onset of negative 
energy free wave solutions (Klein paradox [4,5])? How does particle statistics and in particular the 
Pauli exclusion principle contribute to or modify this. The use of field theory will be mentioned in 
the conclusions of this work. 
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The weak interactions are characterized by the exchange of very massive particles (M^^, Z" ex- 
change) , parity violation that constitute a complicated (but calculable) mix of attractive and repulsive 
forces. For simplicity, one can consider only the Z" exchange. Z° exchange allows for a potential treat- 
ment since it is well represented by a single Feynman graph in momentum space from which a Yukawa 
potential can be derived by Fourier transform. As a consequence of the V-A nature of weak interac- 
tions, it can be shown that according to the total spin state of the two fcrmion system the potential 
is either attractive or repulsive. A similar attractive/repulsive alternative occurs for the isospin state 
in p-n system from which the deuteron singlet bound state emerges. The limitation for the existence 
of an electron-neutrino bound state is set by the enormous value of ^ where ^,{~ 91Gev) is the Z° 
mass [6] and m the reduced neutrino mass. As we have said, oscillation experiments are consistent 
with the existence of neutrino mass states but these have masses of less than a few electron-volts. 

In this paper, we wish to investigate a small part of this question. We concentrate our attention 
upon the Yukawa potential and ask what are the limits upon the exchanged mass for two fcrmion 
bound states to exist. The attractive Coulomb potential yields infinite bound states solutions inde- 
pendent of the coupling strength for both the Schrodingcr and Dirac equations. A Yukawa potential 
will on the contrary not yield a bound state unless the coupling is sufficiently " strong" . This is proven 
only numerically since no analytic solution is known for the Yukawa potential. Using the Schrodinger 
equation and 

y--(r) = -A^, (1) 

(A > for bound states) it is known that the condition for the existence of the lowest lying S-states 

u 

A > 0.84 — . (2) 

m 

An analytic approximation to this result can be found by using the Hulthcn potential [2, 7] in the 
Schrodinger equation. 

The numerator in this potential has been chosen not only because of dimensional requirements but 
also to agree with Eq.JQ) up to the linear term in fir. The imposition of a normalized wave function 
here requires 

\> ^. (4) 
m 

The radial solution for the Hulthcn ground state is 

i?(r, 1 = 0)^2 ^mX{m^X^ ~ fi^) '-^^^^ e"™^'' . 

An equivalent potential (with analytic solution) to the Hulthcn for the Dirac equation is not known. 
We must qualify this last statement. If one allows for vector as well as a scalar potential, and one 
makes a very particular choice then one can find analytic solutions with a Hulthen scalar potential [8] . 
This technique derives from some ingenious suggestions by Alhaidari [9] for solving the Dirac-Morse 
problem. 



II. DIRAC EQUATION AND YUKAWA POTENTIALS 

The Dirac equation in the presence of a general spherical potential V(r) can be written as 



By using 



*(r) 



./i^M3^j™,(r) 



where 3^j„i^. {i") arc the spherical functions obtained by summing the spherical harmonics Y"; mj±^ with 
the spinor eigenstates and k = ±( j + i ) for Z = j±i, for details see ref. [10], Eq.© reduces (dropping 
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the subscripts and superscripts for the radial functions and g^) to two coupled first order ordinary 
differential equations 



[E-m~V{T)] flr) = -g'{r)- 



[E + m-V{r)] 9{r) = J'{r) + 



1 - k 



9ir) , 



(6) 



fir) 



It is instructive to make the non-relativistic limit in Eqs.® by setting E = m— e{e>0 for bound 
states) and assuming e <C m. Eliminating g (the "small" components) yields a second order equation 
for / (we drop both e and V{r) with respect to 2m in the second of Eqs.lO) 



or equivalently 



[e + V{r)] f{r) = 

f'jr) f'jr) 
2m mr 



f'jr) ^ f'jr) fc(fc + l) 
2m mr 2mr^ 

1(1 + 1) 



fir) 



2mr^ 



Vir) 



/(r) = -e/(r) 



This is just the radial part of the Schrodinger equation 



V" 
2m 



Vir) 



(r) = £;„„0(r) 



(7) 

(8) 

(9) 



where Ej,j, = -e < and <p{r) fir)Yirni (r). 

A more sophisticated limit exists, where relativistic corrections are maintained (up to order p"). 
This equation can be derived either by a Foldy-Wouthuysen transformation [11], or in the more 
heuristic method used by Sakurai [12] in which the / function is corrected in order to be normalized 
to order p'^. This equation reads (always assuming a spherical potential so that WV{r) = V'{r) r) 



with 



PL 

2m 



Vir) 



1 



8m''' 8m^ 



Enr ' 



V'V{r) 



1 



Am^r 



V'{r)(T-L 



Now consider the Yukawa potential ^^'"''(y-). For the case I = (S-wave) Ea. (|ll() reduces to 



H 



[i = n] 



P* 



X 



1 



1 

8m^ 



V 



-fir 



2m 8m^ 

The last term in the above equation contains the Darwin delta term 



,2 e 



r 



r 



■ATT5\r) 



(10) 
(11) 

(12) 
(13) 



This term is essential in the Coulomb case (jjL = 0) to reproduce to order p" the bound state energy 
dependence on n (principal quantum number) and j only. When ^ ^ 0, we note that this term contains 
an additional piece proportional to /i^ and hence with the same sign as the original potential term. It 
produces an effective potential with the coupling constant "amplified" 



A 1 



(14) 



If this result is combined with the numerical Schrodinger calculations mentioned previously we might 
expect the condition for the existence of a bound state to be 



(15) 
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If true not only would there be bound states for a given A with values of ^ higher than otherwise 
expected, but far more spectacularly, for ^ ^ 1 the condition for the existence of a bound state is 
"inverted" and reads 

777 

A > 6.62 - (^<1)- (16) 

Of course such a conclusion is highly speculative since it is based upon the combination of Schrodinger 
results and only a part of the relativistic correction to the Schrodinger equation. Higher order cor- 
rections could greatly modify this hypothesis. At this point it is logical to go directly to the fully 
relativistic Dirac equation and (numerically) solve it for the Yukawa potential. We wish to see if there 
is an amplification effect and its comparison with Eg. ljlSI) . 

All our results arc conveniently expressed in terms of the dimensionless parameters 



A, w = ^ and r;=(|-l)/(Vr^-l) 

TO A \ m / 



Using these variables Eqs.© can be written as 



v'{x) 



k 



v{x) 



u'{x) + — u{x) 



\/l-A^-l 



+ A- 



u{x) = 



v{x) = , 



(17) 



where x ~ mXr, u ^ f/x and v — g/x. Note that here A appears explicitly. The Schrodinger limit 
(A <^ 1) is particular in that it can be expressed purely in terms of rj and w with no explicit dependence 
upon A, 



1(1 + 1) e~^^ 
"sch(-^) 1 — Usch(a;) - i]Us,i-,{x) + 2 Msch(a;) =0 



[usch(a;) = f{x)/x] 



(18) 



The asymptotic behavior for x —^ oo and .t — > are substantially different in the Schrodinger and 
Dirac equations. In the former, we have 



Usoh(a; oo) 



while 



Unir{x oo) 



and WDir(a; — ^ oo) 



The Dirac limits constrain the possible values of the bound state energy hy E < m? ov m> E > —m, 
otherwise the solution will not be localized. The lower limit is an example of the Klein paradox, here 
involving a bound state fermion. Indeed, since V{oo) = for E < —to we would have at infinity 
E — V < —TO, which for a step potential is exactly the condition for the Klein paradox . 

There is also a substantial difference between the Schrodinger and Dirac equations in the limit of 
a; — > (recall that x is proportional to the radial parameter). For Schrodinger one has 



X 



while 



Unir{x 0) 



Usoh(a; 0) 
"■^^ and also v^^Jx 



X 



0) 



with = fc^ — A^. Since v must be real and greater than — i for a normalizable solution this sets a 
limit upon the size of A whose value depends upon k i.e. angular momentum. Given that the minimum 
value of is zero, we have that 

A' < F . (19) 



For the S-wave (fc 
exists), we obtain 



-1) which we expect to be the lowest energy bound state (if a bound state 
AlTx°™ < 1 • (20) 
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We also recall the well-known result [10] that since 

V = \J\ - A' < 1 (fc = -1) . 

the Dirac S-wave function is infinite at .t — > 0, both for Yukawa and Coulomb potentials. This is in 
starch contrast with the Schrodinger results. 

In Figs. 1 and 2, we show the numerical results for the S-wave ground state. In Fig. 1 the Dirac 
equation result of r\ versus w is plotted. As the ratio of exchanged to bound state mass ^ is increased 
(increasing w) the value of r\ is reduced. This corresponds to the bound state energy increasing towards 
the limit oi E = m beyond which the bound state ceases to exist. For the case displayed with A = 0.01, 
we essentially reproduce the Schrodinger equation curve, which is A independent as long as A ^ 1. 
This curve confirms the well known condition < 1.19 A or w < 1.19 for Schrodinger. On the same 

m _ . . . 

graph, we display the results of the numerical solutions to the Dirac equation for various values of A. 
As A increases the Dirac results yield higher values of w^^^ and hence evidence for an amplification 
effect. In Fig. 2, we plot an interpolation curve of w^^x against A (up to its maximum S-wave value of 1) 
for the Dirac results. The Schrodinger curve is a flat line at 1.19 on this plot. Again the magnification 
effect is evident. As a comparison, we display in the same plot the prediction of the Darwin term, 
which stimulated this analysis. For A < 0.8 the Dirac amplification is lower than for the Darwin term, 
but it exceeds the latter for A ~ 1. At A = 1 the Schrodinger equation predicts as a limit for an S-wave 
bound state 

(-)""=1.19, (21) 

V TO / max 

with Dirac we find 




while the Darwin term yields 




III. CONCLUSIONS 

Within the validity of the Dirac equation (A^ k^), we have confirmed by numerical calculations that 
the effective Yukawa coupling constant is amplified with respect to the non-relativistic Schrodinger 
potential. However, we have not detected the spectacular "turn-over" phenomena as suggested by the 
Darwin term for the Yukawa potential. At the very least our results imply that bound states exist for 
higher mass exchanges than otherwise expected. We have also identified in the analytic studies of Gu, 
Zheng and Xu [8] the presence of an amplification effect. However, the particular choice of potential 
makes this result of doubtful application. 

We are somewhat unsatisfied by the limits upon coupling constant and bound state energy set by 
the asymptotic Dirac equations. To go beyond this, we must necessarily use field theory. This involves 
a numerical analysis of the Bethe-Salpeter equation [13], and such an analysis is indeed planned. 
Furthermore, the limitation set by the Klein effect (E > —to) is very interesting for a bound state 
fermion. Conventionally, this effect is interpreted by invoking pair production [5]. Since, an attractive 
(binding) potential for say a fermion is repulsive for the corresponding anti- fermion we expect, if pair 
production occurs, an anti-fermion flux to leak from the system while the density of trapped fcrmions 
increases. However, the Pauli exclusion principle will eventually block an increase in the number of 
bound fermions. For example, with Yukawa there may be only one bound state level which could 
accommodate at most two spin ^ fermions. This suggests that the Pauli principle could block the 
Klein effect and allow for bound states with {E < —to). 
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